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Abstract 

We have re-analyzed the world data on inclusive polarized DIS, in both NLO 
and LO QCD, including the very precise JLab Hall A neutron data, and have 
studied the role of positivity constraints, and demonstrated their importance, 
in determining the strange and gluon densities. We have shown that higher 
twist corrections are essential in the analysis of the present data on the structure 
function gi. A consistent QCD analysis is achieved, and results for the polarized 
parton densities are given in both the MS and JET schemes. 


1 Introduction 


Spurred on by the famous European Muon Collaboration (EMC) experiment [1] at 
CERN in 1987, there has been a hnge growth of interest in the partonic spin structnre 
of the nncleon, i.e., how the nncleon spin is bnilt np ont from the intrinsic spin and 
orbital angnlar momentnm of its constituents, qnarks and glnons. Onr present knowl¬ 
edge about the spin strnctnre of the nncleon comes from polarized inclusive and semi- 
inclusive DIS experiments at SLAC, CERN, DESY and JLab, polarized proton-proton 
collisions at RHIC and polarized photoprodnction experiments. The determination of 
the longitndinal polarized parton densities in QCD is one of the important aspects of 
this knowledge. 

In this paper we present an npdated version of our NLO QCD polarized parton 
densities in both the MS and the JET (or so-called chirally invariant) [2] factorization 
schemes, as well as the LO ones, determined from the world data [1, 3, 4] on inclusive 
polarized DIS. Comparing to onr previous analyses [5, 6]; i) The recent very precise 
JLab/ Hall A data [4] on Qi /are incorporated into the analysis and ii) New positivity 
constraints are imposed and their role in the determination of the polarized parton 
densities is discussed. The updated polarized parton densities (PPD) are compared 
to those obtained by the other groups and the effect of different positivity constraints 
that have been imposed is demonstrated. 

2 QCD fits to the data 

In QCD the spin strnctnre fnnction gi can be written in the following form (Q^ >> A^): 

gi{x, Q"^) = gi{x, Q^)lt + gi{x., Q^)ht , (1) 

where ”LT” denotes the leading twist (r = 2) contribntion to gi, while ”HT” denotes 
the contribntion to gi arising from QCD operators of higher twist, namely r > 3. In 
(1) we have dropped the nncleon target label N. The HT power corrections (np to 
Oi^l/Q"^) terms) can be divided into two parts: 

gi{x, Q^)ht = h™^{x, Q^)IQ^ F h[x, Q^)IQ^ , (2) 

where are the calculable [7] kinematic target mass corrections and effec¬ 

tively belong to the LT term. h{x, Q'^) are the dynamical higher twist (r = 3 and 
r = 4) corrections to gi, which are related to mnlti-parton correlations in the nucleon. 
The latter are non-pertnrbative effects and cannot be calculated without using models. 
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(Note that the twist-3 contribution to gi is equal to zero if the Wandzura-Wilczek 
approximation [8] for the spin structure function g 2 is used.) (5 ^)lt in (1) is the 

well known pQCD expression and in NLO has the form 


gi { x , Q^)pqcd “ 2 ^ ® (IH— 


5C,) + 




where Ag(a;, Ag(a;, and AG(a;, are quark, anti-quark and gluon polarized 
densities in the proton, which evolve in according to the spin-dependent NLO 
DGLAP equations. 5C{x)q^G are the NLO spin-dependent Wilson coefficient functions 
and the symbol ® denotes the usual convolution in Bjorken x space. Nf is the number of 
active flavors. Note that in the LO QCD approximation the coefficients 5C{x)q ,G in (3) 
vanish and the polarized parton densities evolve according the LO DGLAP equations. 

One of the features of polarized DIS is that a lot of the present data are in the 
preasymptotic region ~ 1 — 5 GeG^,4 GeV'^ < < 10 GeV"^). While in the 

unpolarized case we can cut the low and data in order to minimize the less well 
known higher twist effects, it is impossible to perform such a procedure for the present 
data on the spin-dependent structure functions without losing too much information. 
This is especially the case for the HERMES, SLAG and Jefferson Lab experiments. So, 
to confront correctly the QGD predictions with the experimental data and to determine 
the polarized parton densities special attention must be paid to the non-perturbative 
higher twist (powers in 1/Q^) corrections to the nucleon structure functions. 

We have used two approaches to extract the polarized parton densities from the 
world polarized DIS data. According to the hrst [5] the leading twist LO/NLO QGD 
expressions for the structure functions gi and Fi have been used in order to confront 
the data on spin asymmetry Ai(~ fi'i/Ei) and gi/Fi. We will refer to these as ^gi/Fi 
hts. We have shown [9, 10] that in this case the extracted from the world data ’effective’ 
HT corrections to the ratio gi/Fi 


9i{x,Q^)'\ _ gi{x,Q^)LT 

F,{x,Q^)\^^^ F,{x,Q^)^t 

are negligible and consistent with zero within the errors, i.e. ~ 0, when 

for {gi)LT and (Ei)lt their NLO QGD approximations are used. (Note that in QGD 
the unpolarized structure function Fi takes the same form as gi in (1), namely Fi = 
+ {Fi)ht-) What follows from this result is that the higher twist corrections to 
gi and Fi approximately compensate each other in the ratio gi/Fi and the NLO PPDs 
extracted this way are less sensitive to higher twist effects. This is not true in the LO 
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case (see our discussion in Ref. [6]). The set of polarized parton densities extracted 
this way is referred to as or PD(Set 1). 

According to the second approach [6], the gi/Fi and Ai data have been fitted 
using phenomenological parametrizations of the experimental data for the unpolarized 
structure function ^ 2 ( 2 ^, Q'^) and the ratio R{x, Q^) of the longitudinal to transverse yiV 
cross-sections (i.e. Fi is replaced by its expression in terms of usually extracted from 
unpolarized DIS experiments F 2 and R). Note that such a procedure is equivalent 
to a fit to {gi)exp, but it is more consistent than the fit to the gi data themselves 
actually presented by the experimental groups because here the gi data are extracted 
in the same way for all of the data sets. In this case the HT corrections to gi cannot 
be compensated because the HT corrections to ^ 1(^2 and R) are absorbed in the 
phenomenological parametrizations of the data on F 2 and R. Therefore, to extract 
correctly the polarized parton densities from the gi data, the HT corrections (2) to 
gi have to be taken into account. In our £t to the data we have used the following 
expressions for gi/Fi and Ai: 


Fr^{x,Q^) 

A^{x,Q^) 


exp 


exp 






g^{x, Q^)lt + h^{ x)/Q‘^ [1 -h R{x, Q^)exp] 


F^{x,Q^] 


exp 


^f(a;,Q2)LT + 


Fi^{x,Q‘^] 


■2a;[l -|- R{x, 


I expj ? 


( 5 ) 


exp 


where gi{x, Q^)lt (N=p, n, d) is given by the leading twist expression (3) in LO/NLO 
approximation including the target mass corrections through the Nachtmann variable 
[11, 12]. In (5) h^{x) are the dynamical r = 3 and r = 4 HT corrections which are 
extracted in a model independent way. In our analysis their dependence is neglected. 
It is small and the accuracy of the present data does not allow to determine it. For the 
unpolarized structure functions F^{x,Q'^)exp and R{x,Q^)exp we have used the NMC 
parametrization [13] and the SLAG parametrization Riggs [14], respectively. We will 
refer to these as \gi + HT)’ fits, and the set of polarized parton densities extracted 
according to this approach - -|- HT) or PD(Set 2). 

As in our previous analyses [5, 6], for the input LO and NLO polarized parton 
densities a.t Ql = 1 GeV^ we have adopted a simple parametrization 


xAuy{x,Ql) = 
xAd^{x,Ql) = 
xAs{x,Ql) = 
xAG{x,Ql) = 


r]uAuX°‘'‘xUy{x,Ql), 

ridAdx''‘^xdy{x,Ql), 

risAsX''‘’xs{x,Ql), 

r]gAgx‘"<’xG{x,Ql), 


( 6 ) 
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where on RHS of (6) we have used the MRST98 (central gluon) [15] and MRST99 
(central gluon) [16] parametrizations for the LO and NLO(MS) unpolarized densities, 
respectively. The normalization factors Ai in (6) are fixed such that rji are the first 
moments of the polarized densities. To £t better the data in LO QCD, an additional 
factor (l+ 7 ^a;) on the RHS is used for the valence quarks. Bearing in mind that the light 
quark sea densities Ah and Ad cannot, in principle, be determined from the present 
inclusive data (in the absence of polararized charged current neutrino experiments) we 
have adopted the convention of a flavor symmetric sea 

AUsea = Au = Adsea = Ad = As = As. (7) 

The first moments of the valence quark densities rju and 77 ^ are constrained by the 
sum rules 

ag = = F + D = 1.2670 ± 0.0035 [17], (8) 

as = 3F - D = 0.585 ± 0.025 [18], (9) 

where Og and Og are non-singlet combinations of the first moments of the polarized 
parton densities corresponding to 3“'^ and 8**^ components of the axial vector Cabibbo 
current 

as = (Au + Ah)(g2) _ (Ad + Ad) (Q^) , (10) 

as = (Au + Ah)(g^) -h (Ad -h Ad){Q‘^) - 2{As + As){Q‘^) . (11) 

The polarized parton densities (6) and (7) have to satisfy the positivity condition, 
which in LO QCD implies: 

\Afi{x,Ql)\ < \Afi{x,Ql)\ < fi{x,Ql). (12) 

The constraints (12) are the consequence of a probabilistic interpretation of the 
parton densities in the naive parton model, which is still valid in LO QCD. Beyond 
LO the parton densities are not physical quantities and the positivity constraints on 
the polarized parton densities are more complicated. They follow from the positivity 
condition for the polarized lepton-hadron cross-sections Adj in terms of the unpolarized 
ones (|A(jj| < aQ and include also the Wilson coefficient functions. It was shown [19], 
however, that for all practical purposes it is enough, at the present stage, to consider LO 
positivity bounds for LO as well as for for NLO parton densities, since NLO corrections 
are only relevant at the level of accuracy of a few percent. 

While in our previous NLO QCD analyses we have mainly used for the unpolarized 
parton densities on the RHS of (12) the Barone et al. parametrization [20] we are here 
using the MRST02(NLO) updated unpolarized parton densities [21] in both MS and 
JET schemes. The only significant change in the MRST02 NLO partons, compared to 
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Table 1. The parameters of the Set 1 of NLO input parton polarized densities 
at = 1 as obtained from the hts to the world [1, 3] and 

JLab [4] data in the MS and JET schemes. The errors shown are total (statistical and 
systematic). The parameters marked by (*) are hxed by the sum rules (8) and (9). 


Fit 



DF 

188 - 6 

188 - 6 


159.0 

158.5 

XVDF 

0.874 

0.871 

Vu 

0.926* 

0.926* 

(Xu 

0.209 ± 0.018 

0.211 ± 0.018 

Vd 

- 0.341* 

-0.341* 


0.072 ± 0.068 

0.083 ± 0.058 

Vs 

- 0.066 ± 0.009 

- 0.049 ± 0.013 

(Is 

0.649 ± 0.117 

0.779 ± 0.214 

Vg 

0.195 ± 0.257 

0.289 ± 0.316 

ttg 

2.575 ± 1.729 

0.000 ± 0.686 


those of MRST99 (which we have also used in some of our analyses), is in the gluon, 
and especially in an increase in the gluon density at high x. In the MRST02 fit to 
the world data the authors no longer include prompt photon data due to theoretical 
problems and possible inconsistencies between data sets, and instead allow the high 
X gluon to be determined by the improved Tevatron jet data [22], which considerably 
improves the determination of the gluon. The use of new positivity constraints leads to 
a signihcant change in the polarized strange quark and gluon densities, the importance 
of which will be discussed further. 


3 Results 

In this section we present the numerical results of our hts to the world data [1, 3] on 
gi/Fi and Ai including JLab Hall A neutron data [4]. The data used (188 experimental 
points) cover the following kinematic region: 

0.005 < a; < 0.75, 1 < < 58 GeV^ . (13) 
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Table 2. The parameters of the Set 2 of LO, NLO(MS) and NLO(JET) input 
parton densities PD( 5 fj"'^ + HT) at = 1 GeV^ as obtained from the best {gi + HT) 
fits to the world [1, 3] and JLab [4] data. The errors shown are total (statistical and 
systematic). The parameters marked by (*) are hxed. Note that the TMC are included 
in i9i)LT- 


Fit 

(fi-OLO + Kx)/Q‘^ 

(5'i)nlo(ms) + h{x)/Q‘^ 

(fl'l)NLO(JET) + h{x)/Q‘^ 

DF 

188 - 16 

188 - 16 

188 - 16 


152.7 

150.0 

150.8 

XVDF 

0.888 

0.872 

0.877 

Vu 

0.926* 

0.926* 

0.926* 

CLu 

0.000 ± 0.011 

0.244 ± 0.036 

0.240 ± 0.037 

In 

1.607 ± 0.292 

0* 

0* 

Vd 

- 0.341* 

-0.341* 

-0.341* 

Q-d 

0.000 ± 0.067 

0.123 ± 0.134 

0.120 ± 0.148 

Id 

3.184 ± 1.630 

0* 

0* 

Vs 

- 0.075 ± 0.010 

- 0.078 ± 0.012 

- 0.064 ± 0.017 

Q/g 

0.514 ± 0.068 

0.629 ± 0.090 

0.612 ± 0.121 

Vg 

0.602* 

0.348 =h 0.345 

0.268 ± 0.422 

ttg 

0.328* 

1.980 ± 1.359 

2.851 ± 1.444 

Xi 

hP{xi) [GeV^] 

0.028 

-0.003 ± 0.036 

0.010 ± 0.042 

0.024 ± 0.039 

0.100 

- 0.097 ± 0.033 

- 0.043 ± 0.034 

- 0.040 ± 0.039 

0.200 

- 0.164 ± 0.032 

- 0.106 ± 0.036 

- 0.108 ± 0.038 

0.350 

- 0.036 ± 0.036 

- 0.016 ± 0.038 

- 0.015 ± 0.040 

0.600 

0.032 ± 0.019 

0.046 ± 0.019 

0.048 ± 0.019 

Xi 

h^{xi) [GeV^] 

0.028 

0.204 ± 0.078 

0.145 ± 0.081 

0.161 ± 0.081 

0.100 

0.168 ± 0.050 

0.192 ± 0.047 

0.197 ± 0.046 

0.200 

0.023 ± 0.058 

0.035 ± 0.067 

0.032 ± 0.070 

0.325 

0.031 ± 0.027 

0.019 ± 0.031 

0.018 ± 0.035 

0.500 

0.031 ± 0.015 

0.010 ± 0.014 

0.011 ± 0.016 
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The total (statistical and systematic) errors are taken into acconnt. The systematic 
errors are added qnadratically. 

We have determined from the data two sets of polarized parton densities in both the 
MS and the JET factorization schemes: (see Table 1) and PD( 5 f['"'" + 

HT) (see Table 2). For the second set of polarized parton densities a LO version 
PD( 5 f^*^ + HT) is also presented. The latter have been extracted from the data nsing 
for the nnpolarized parton densities on the RHS of the positivity bonnds (12) the LO 
MRST’Ol ones (the MRST’02 LO set is essentially identical to the MRST’Ol one and 
is not presented by the anthors). Note once more that in the LO ’’gi/Fi fit to the data 
the polarized parton densities PD( 5 f}"®/F^'"®) cannot be correctly determined becanse 
the Callan-Gross relation 2a;Fi(a;, (5 ^)lo = ^ 2 ( 2 ^, Q^)lo (used in the calcnlation of Fi) 
is strongly broken (by np to 30 %) for x < 0.25 and small Q^. In other words, to 
extract correctly the LO polarized parton densities from the data, the second method 
of analysis have to be nsed, i.e. the high twist corrections to gi have to be taken into 
acconnt. This observation is especially important for the analysis of the semi-inclnsive 
DIS data, where the LO QCD approximation is mainly used. 




Figure 1: NLO(MS) polarized parton densities PD( 5 fj""'" + HT) (solid curves) together 
with their error bands compared to PD( 5 fJ^^®/Fj^^°) (dashed curves) at = 4 GeV^. 

In Fig. 1 we compare the NLO(MS) polarized parton densities PD( 5 f{'"'" + HT) with 
PD( 5 fJ^^®/Fj^^*^). As seen from Fig. 1 the two sets of polarized parton densities are 
very close to each other, especially for u and d quarks. This is a good illustration of 
the fact that a fit to the gi data taking into account the higher twist corrections to gi 
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(Xdfnlo = 0.872) is equivalent to a fit of the data on Ai(~ gi/Fi ) and gi/Fi using 
for the gi and Fi structure functions their NLO leading twist expressions (xdf nlo = 
0.874). In other words, this fact conhrms once more that the higher twist corrections 
to gi and Fi approximately cancel in the ratio gi/Fi. Nevertheless, we consider that 
the Set 2 of the polarized parton densities + HT) is preferable because using 

them and simultaneously extracted higher twist corrections to gi, the spin structure 
function gi can be correctly calculated in the preasymptotic (Q^, hh^) region too. 




Figure 2: Comparison between NLO Set 2 polarized strange sea and gluon densities 
at = 4 GeV^ in the MS (solid curves) and the JET scheme (dashed curves). Dot 
curve corresponds to As(a;, Q^)jet^ms obtained by the transformation rule (14) (see 
the text). 

Let us briefly discuss the scheme dependence of our results. Recall that according 
to perturbative QCD the NLO polarized valence quarks and gluons should be the same 
[up to order of 0{al)] in both the MS and the JET factorization schemes, while for 
the moments of the strange sea quarks the following transformation rule is valid: 

As(n, Q^)ms = As(n, Q^)jet - 2Tin(^n + 1) • (14) 

It is seen from the Table 1 and Table 2 that the values of x^/DF(MS) and 
X^/DF(JET) coincide almost exactly for the gi/Fi as well as for the {gY^ + HT) hts, 
which is a good indication of the stability of the analysis regardless of the scheme used. 
To illustrate the factorization scheme dependence, the extracted polarized PD( 5 ff^*^ + 
HT) in the schemes under consideration are compared in Fig. 2. Note that the valence 
densities An^ and Ad^ in the MS and JET schemes are almost identical (see the 
values of the parameters for the corresponding input parton densities in Table 2) and 
in excellent agreement with what follows from QCD. So the corresponding curves are 
not shown in Fig. 2. The extracted polarized gluons in the two schemes are also 





well consistent within the errors (see Fig. 2). In Fig. 2 we also show the polarized 
strange sea densities in both schemes, determined directly from the hts and evolntion 
eqnations, as well as the strange sea As{x, Q^)jet^ms obtained by the transformation 
rule (14). It is seen that: i) at large x As(a;, (5 ^)ms is very close to As(a;, (5 ^)jet 
which is consistent with (14), since large n in the Mellin space corresponds to large 
Bjorken x, and ii) As(a;, <5 ^)jet=^ms coincides very well with As(a;, <5 ^)ms. We have 
discussed here in detail the results on the set PD( 5 f^^° + HT). We have found the 
same conclusion for the other set of polarized parton densities PD( 5 fJ^'"®/F)^^*^). In 
conclusion, we have found that the obtained numerical results are in a good agreement 
with the perturbative QCD predictions. 

Compared to our previous results [5, 6] we now obtain smaller values for the gluon 
polarization (the hrst moment of AG{x,Q‘^)), which leads to a smaller difference be¬ 
tween the values of the strange quark polarization (the hrst moment of As{x,Q^)) 
determined in the MS and the JET schemes, respectively. To illustrate this tendency 
we present in Table 3 the corresponding values of the hrst moments. As, AG and AS, 
for the LSS’Ol and LSS’05(Set 1) sets of polarized parton densities. Note that in the 
JET scheme the singlet polarization AS((5^) is a independent quantity. Then, in 
this scheme it is meaningful to directly interpret AS as the contribution of the quark 
spins to the nucleon spin and to compare its value obtained in the DIS region with the 
predictions of the diherent (constituent, chiral, etc.) quark models at low Q^{Q^ ~ A^). 
Our new value of ASjet = 0.29 ± 0.08 is smaller then the old one and farther from 
the value 0.6 of AS at low region predicted in relativistic constituent quark models 
[23]. Note, however, that if nonperturbative vacuum spin ehects are taken into account 
[24, 25], the value of AS at low is expected to be smaller than 0.6. Further theoreti¬ 
cal and experimental investigation in both the large and very low regions, is needed 
to answer more precisely the important question as to the fraction of the nucleon spin 
carried by its quarks. 


Table 3. First moments (polarizations) of LSS’Ol and LSS’05(Set 1) polarized 
parton densities at = 1 GeV^. 


Fit 


^'S(Q^)jET 

AG((5^)jet 

ASjet 

LSS’Ol 

-0.065 ± 0.016 

-0.035 ± 0.010 

0.68 ± 0.32 

0.37 ± 0.07 

LSS’05/Set 1 

-0.066 ± 0.009 

-0.049 ± 0.013 

0.29 ± 0.32 

0.29 ± 0.08 
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The extracted higher twist corrections to the proton and nentron spin structure func¬ 
tions, h^{x) and h^{x), are shown in Fig. 3.f As seen from Fig. 3 the size of the 




Figure 3; Higher twist corrections to the proton and neutron gi structure functions 
extracted from the data on gi in NLO(MS) QCD approximation for (5 ^)lt- The 
parametrization (15) of the higher twist values is also shown. 

HT corrections is not negligible and their shape depends on the target. In Fig. 3 our 
previous results on the higher twist corrections to gi (before the JLab Hall A data were 
available) are also presented. As seen from Fig. 3, thanks to the very precise JLab Hall 
A data at large x the higher twist corrections to the neutron spin structure function 
are now much better determined in this region. In Fig. 3 our parametrizations of the 
values of higher twists for the proton and neutron targets 


h^{x) 

0 1913 

- 0.0465 _ exp[ 2{{x 0.2087)/0.2122)2] 

A/2 


h^{x) 

0 2420 

- 0.0119+ _ exp[ 2{{x 0.0783)/0.1186)2] 

A/2 

(15) 


are also shown. These should be helpful in a calculation of the nucleon structure 
function gi for any x and moderate in the experimental region, where the higher 
twist corrections are not negligible. The values of the higher twist corrections to the 
proton and neutron gi structure functions extracted in a model independent way from 
polarized DIS data are in agreement with the QCD sum rule estimates [27] as well as 
with the instanton model predictions [28] but disagree with the renormalon calculations 

[ 29 ]^_ 

tThe moments of higher twist contribution to the proton structure function gi, including also in 
the analysis the data on the resonance region, have been determined in the very recent paper [26]. 
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Figure 4: Effect of the COMPASS data on the Set 2 NLO(MS) polarized parton den¬ 
sities + HT). 

When this analysis was hnished, the COMPASS Collaboration at CERN reported 
new data on the longitudinal asymmetry Af [30]. Their results improve considerably 
the statistical accuracy on the small x region 0.004 < x < 0.03. We have carried out 
a preliminary study {'gi -|- HT’ fit in the MS scheme) to see whether the COMPASS 
data has any signihcant effect on the result of our analysis. The effect of the new data 


NLO MS 

^ 0-3 

a 

S 0.1 
0.0 
- 0.1 
0.3 
0.2 
0.1 
0.0 
- 0.1 

0.0 0.2 0.4 0.6 0.8 

X 

Figure 5: Effect of the COMPASS data on the higher twist values. 

on the polarized parton densities PD( 5 ff^*^ -|- HT) and the higher twist corrections is 
illustrated in Fig. 4 and Fig. 5. While the valence quarks densities and do 
not change in the experimental region (for that reason they are not shown in Fig. 4), 
the magnitudes of both the polarized gluon and strange quark sea densities decrease, 
but the corresponding curves he within the error bands (see Fig. 4). The impact of the 
new data on the values of higher twist corrections is negligible (see Fig. 5). The new 


• World data+ JLab/HallA/n 

O World data+ JLab/HallA/ti - 
+ COMPASS/d 

1 ■ 


I" Proton. 



[ 

Neutron- 
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values are in good agreement with the old ones although there is a tendency for the 
central values for the proton target to be slightly lower than the old ones. The central 
values of the HT correction for the neutron at small x are also slightly lower than the 
old ones. 


4 Impact of positivity constraints on polarized PD 

Let us consider now how the use of different positivity constraints influences the results 
on the polarized parton densities. In Fig. 6 we compare our new Set 1 of NLO(MS) po¬ 
larized parton densities with LSS’2001 parton densities [5] presented 

on the HEPDATA web site. Both sets are determined from the data by the same 
method but using different positivity constraints. (Note that the inclusion of the JLab 
data do not influence the results on the polarized parton densities if the same positivity 
constraints are used.) While the new polarized are compatible with 

the positivity bounds (12) imposed by the MRST’02 unpolarized parton densities [21], 
those of the LSS’2001 set are limited by the Barone et ah unpolarized parton densities 
[20]. As seen from Fig. 6 the valence quark densities and Ad^ of the two sets are 
close to each other, while the polarized strange sea quark and gluon densities are sig- 
nihcantly different. This comparison is a good illustration of the fact that the present 
inclusive polarized DIS data allow a much better determination of the valence quark 




Figure 6: Comparison between our two sets of NLO(MS) polarized parton densities, 
LSS’Ol and LSS’05(Set 1), at = 4 GeCT 
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X 


Figure 7: Comparison between the NLO(MS) unpolarized strange quark sea and gluon 
densities determined by MRST’02 [21] and Barone at al. [20]. 

densities (if SU(3) symmetry of the flavour decomposition of the sea is assumed) than 
the polarized strange quarks As{x,Q‘^) and the polarized gluons AG{x,Q^). This is 
especially true for the high x region, where the values of As{x, Q"^) and AG{x, Q^) are 
very small and the precision of the data is not enough to extract them correctly. That 
is why different unpolarized sea quark and gluon densities (see Fig. 7) used on the RHS 
of the positivity constraints (12) are important and crucial in determining As(x, Q^) 






Figure 8: Comparison between our NLO(MS) polarized parton densities (Set 1) at 
= 4 GeV'^ with those obtained by GRSV (’standard scenario’) [31], BB (ISET=4 
or BB2) [32] and AAC (AAC03) [33]. 

and AG{x, Q^) in this region. The more restrictive s(x, Q^)mrst '02 at high x leads to a 
smaller value of |As(x, Q^)|lss '05 ia this region, while the smaller G(x, Q^)Bar.et.ai pro- 
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vides a stronger constraint on AG(a;, (5^)lss'oi (see Fig. 6). To illustrate this fact once 
more, we compare our new polarized parton densities at = 4 GeV^^ 

with those obtained by GRSV [31], Blumlein, Bottcher [32] and the Asymmetry Analy- 




Figure 9: Comparison between the NLO(MS) unpolarized strange quark sea and gluon 
densities determined by MRST [16, 21] and GRV [34], 

sis Collaboration (AAC) [33] using almost the same set of data (see Fig. 8). Note that 
all these groups have used the GRV unpolarized parton densities [34] for constraining 
their polarized parton densities at large x. As seen from Fig. 9, in this x region the 
unpolarized GRV and MRST’02 gluons are practically the same, while the magnitude 
of the unpolarized GRV strange sea quarks is much smaller than that of MRST’02. 
Therefore, the GRV unpolarized strange sea quarks provide a stronger constraint on 
the polarized ones. The impact on the determination of the polarized strange sea 
density is demonstrated in Fig. 10. (The GRV and MRST’02 unpolarized strange sea 
densities are also shown.) As a result the magnitude of our polarized strange sea density 
x\As{x,Q‘^)\ is larger in the region x > 0.1 than those obtained by the other groups. 
Note also that the magnitude of xAs obtained by the GRSV and BB is smaller than 
that determined by AAC. We consider the GRSV result to be a consequence of the fact 
that in their analysis, the GRV positivity constraint is imposed at lower value of Q^: 

= j^NLO = 0.4 GeV^, while AAC has used the same requirement at = 1 GeV^. 
Finally, the different positivity conditions on As influence also the determination of the 
polarized gluon density for larger because the evolution in mixes the polarized 
sea quarks and gluons. To end this section we would like to emphasize that for the 
adequate determination of polarized strange quarks and gluons at large x, the role of 
the corresponding unpolarized densities is very important. That is why the latter have 
to be determined with good accuracy at large x in the preasymptotic (Q^, W^) region 
too. Usually the sets of unpolarized parton densities, presented in the literature, are 
extracted from the data on DIS using cuts in and chosen in order to minimize 
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the higher twist effects. In order to use the densities for constraining the polarized 
parton densities they have to be continued to the preasymptotic W^) region. It is 
not obvious that the continued unpolarized parton densities would coincide well with 
those obtained from the data in the region (Q^ > 1 GeV^, > 4 GeV^) in the 
presence of the HT corrections to unpolarized structure functions Fi and F 2 . So, a 
QCD analysis of the unpolarized world data including the preasymptotic (Q^, W^) 
region and taking into account HT corrections is needed in order to extract correctly 
the unpolarized parton densities in the preasymptotic region. Our arguments for the 



Figure 10: Comparison between our NLO(MS) polarized strange sea quark density (Set 
1) at = 1.25 GeV^ with those obtained by GRSV (’standard scenario’) [31], BB 
(ISET=4 or BB2) [32] and AAC (AAC03) [33]. The unpolarized MRST02 and GRV98 
strange sea quark densities are also shown. 

need for a precise determination of the unpolarized densities of strange quarks and 
gluons in both the asymptotic and preasymptotic regions in and W^, coming from 
spin physics, could be considered as additional to those discussed in the recent paper 
|35]. 

5 Conclusions 

i) We have re-analyzed the world data on inclusive polarized deep inelastic lepton- 
nucleon scattering in leading and next-to-leading order of QGD, adding to the old set 
of data the very precise JLab Hall A neutron data. Gompared to our previous analyses 
new positivity constraints on the polarized parton densities have been used. The latter 
reflect mainly the better determination of unpolarized gluon density, especially at high 
X. Two new sets of NLO polarized parton densities in the JET and MS factorization 
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schemes as well as LO polarized patron densities have been extracted from the data 
using different methods of analysis. The NLO polarized parton densities determined 
in the two schemes are in a good agreement with the pQCD predictions. 

ii) The impact of positivity constraints on the polarized parton densities has been 
studied. Special attention has been paid to the role of positivity constraints in deter¬ 
mining the polarized strange quark and gluon densities, which are not well determined 
from the present inclusive DIS data. For that reason the effect of the positivity con¬ 
ditions used to constrain them is very important. On other hand, the different sets of 
unpolarized parton densities needed to impose the positivity bounds (12) are usually 
determined in the DIS region and their continuation to the preasymptotic region via 
the evolution DGLAP equations sometimes leads to very different behaviour of the 
unpolarized parton densities belonging to the different sets. In particular, it is demon¬ 
strated that the use of MRST’02 and GRV98 NLO (MS ) unpolarized strange quark 
densities in the RHS of the positivity condition (12) leads to a signihcant difference of 
the extracted polarized strange sea density, especially for x > 0.05. So, a more precise 
determination of the unpolarized parton densities in the preasymptotic region is very 
important for a better determination of the polarized ones, especially of the polarized 
strange quark sea and gluon densities, which are weakly constrained by the present 
experimental data. 

Hi) It was demonstrated that in the £t to the gi data the higher twist corrections 
to gi are important and have to be taken into account. It was also shown that the 
values of higher twist corrections to the neutron spin structure function at high x are 
determined much more precisely when the JLab Hall A data are used in the analysis. 

iv) Finally, the use of the very recent GOMPASS data on inclusive asymmetry 
Af has little or no effect on the polarized valence densities, but the magnitudes of the 
strange quark and gluon densities decrease. 
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Appendix 

For practical purposes we present here explicitly our Set 1 and Set 2 of polarized 
parton densities at = 1 GeV^. The polarized valence quark densities correspond to 
SU(3) flavour symmetric sea. 

LSS’05 (Set 1) - NLO(MS) PD(c/i/Fi): 


xAuy{x) 

= 0.4621 (1 - ( 

1 + 2.179 a;^/^ + 14.57 a; ) , 


xAdv{x) 

= -0.02257 (1 - 

( 1 + 35.47 a;^/^ + 28.97 a; ) ^ 

) 

xAs{x) 

= -0.02520 (1 - 

( 1 + 3.656 a;^/^ + 19.50 a; ) ^ 

) 

xAG{x) 

= 532.3 (1 - a;)®'®^® ( 1 

-3.147 a;^/2 + 3.148 a; ) . 

(A.l) 


LSS’05 (Set 1) - NLO(JET) PD(c/i/Fi): 


xAuv{x) 

= 0.4635 (1 - ( 

1 + 2.179 a;^/2 + 14.57 a; ) 

xAdv{x) 

= -0.02322 (1 - 

( 1 + 35.47 a;^/2 + 28.97 a; 

xAs{x) 

= -0.02867 (1 - 

( 1 + 3.656 x^/^ + 19.50 a; 

a;AG(a;) 

= 5.080 (1 - a;)®'®^^ ( 1 

-3.147 a;^/2 +3.148 a; ) . 


LSS’05 (Set 2) - LO PD(c/i + HT); 

xAu^{x) = 0.1761 (1 - (1 + 1.607 a;)( 1 - 0.4085 x^^^ + 17.60 a; ) , 

xAd^{x) = -0.00807 (1 - (1 + 3.184 a;)( 1 + 37.25 x^/^ + 31.14 a; ) , 

a;As(a;) = -0.04464 (1 - ( 1 - 0.9052 x^/^ + 11.53 a; ) , 

a;AG(a;) = 1.164 (1 - ( 1 - 4.255 a;^/^ + 7.274 a; ) . (A.3) 

LSS’05 (Set 2) - NLO(MS) PD(^i + HT); 

a;An^(a;) = 0.4958 (1 - ( 1 + 2.179 x^^^ + 14.57 a; ) , 

a;Ad„(a;) = -0.02567 (1 - ( 1 + 35.47 x^^^ + 28.97 a; ) , 

a;As(a;) = -0.02756 ('i _ ^y.(i49 ^ ^ g ggg ^ 1/2 ^ ^ ^ ^ 

a;AG(a;) = 421.9 (1 - ( 1 - 3.147 x^^^ + 3.148 a; ) . (A.4) 
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LSS’05 (Set 2) - NLO(JET) PD(^i + HT); 

xAu^{x) = 0.4924 (1 - ( 1 + 2.179 x^/"^ + 14.57 a; ) , 

xM^{x) = -0.02549 (1 - ( 1 + 35.47 a;^/^ + 28.97 a; ) , 

a;As(a;) = -0.02115 (1 - ( 1 + 3.656 a;^/^ + 19.50 a; ) , 

a;AG(a;) = 1009.5 a;^-®^° (1 - a;)®'®^^ ( 1 - 3.147 a;^/^ + 3.148 a; ) . (A.5) 
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